One fewer solution to the cosmological lithium problem 
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Data from a recent 9 Be( 3 He,£) 9 B measurement are used to rule out a possible solution to the 
cosmological lithium problem based on conventional nuclear physics. 
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The primordial abundance of 7 Li inferred from obser- 
vational data is roughly a factor of three below the abun- 
dance predicted by the standard theory of big bang nu- 
cleosynthesis (BBN) plusing the baryon-to-photon ratio 
rj = 6.19(15) x 10 -10 [2J determined mainly from mea- 
surements of the cosmic microwave background radia- 
tion. In contrast, there is good agreement for 2 H and 
4 He. Taking into account the estimated uncertainties 
on the observationally inferred and the theoretically de- 
duced 7 Li abundances, the significance of the discrepancy 
is 4.2o--5.3cr Q. This constitutes one of the important 
unresolved problems of present-day astrophysics and is 
termed the cosmological lithium problem. Among other 
possibilities, the discrepancy could be due to new physics 
beyond the Standard Model of particle physics [3| , errors 
in the observationally inferred primordial lithium abun- 
dance 1 , or incomplete nuclear physics input for the BBN 
calculations. The present paper addresses the last possi- 
bility. 

In standard BBN theory, assuming 77 = 6.19(15) x 
10~~ 10 , most 7 Li is produced in the form of 7 Bc. 
Only much later, when the Universe has cooled suffi- 
ciently for nuclei and electrons to combine into atoms, 
does 7 Be decay to 7 Li through electron capture. The 
temperature range of 7 Be production is T ~ 0.3 — 
0.6 GK, where the main mechanism for 7 Be production is 
3 Hc(a, 7) 7 Bc while the main mechanism for 7 Be destruc- 
tion is 7 Be(n,p) 7 Li followed by 7 Li(p, a) 4 Hc. The rates 
of these reactions as well as the reactions that control 
the supply of neutrons, protons, 3 He, and a particles are 
known with better than 10% precision at BBN tempera- 
tures Q, resulting in an uncertainty of only 13% on the 
calculated 7 Li abundance [|[. 

A recent theoretical paper Q explores the possibility 
of enhancing 7 Be destruction through resonant reactions 
withp, d, t, 3 He, a, leading to compound states in 8 B, 9 B, 
10 B, 10 C, 11 C, respectively. The paper concludes that, of 
the known excited states in these isotopes [l(| , only 
the 16.8 MeV state in 9 B has the potential to significantly 
influence 7 Be destruction 2 . (Note that in Ref. Q this 
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1 Lithium may be destroyed in metal-poor stars through diffusion 
and turbulent mixing [3]. 

2 Ref. § offers a more optimistic view, but only by adopting a 
somewhat flexible approach to basic principles of nuclear physics. 
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FIG. 1: Schematic illustration of the proposed 7 Be destruc- 
tion mechanism, d + 7 Be — > 9 B* — > p + 8 Be*. The energies 
are in MeV relative to the ground state of 9 B. Subsequently, 
8 Be* breaks up into two a particles. 



state is referred to as the 16.7 MeV state.) The proposed 
destruction mechanism is shown schematically in Fig. [T] 
The 16.8 MeV state in 9 B is formed by the fusion of 
7 Be with a deuteron and decays by proton emission to a 
highly excited state in 8 Be, 16.626 MeV above the ground 
state, which subsequently breaks up into two a particles. 
(The last step is not shown in the figure.) The reason 
why the decay must proceed by proton emission to the 
16.626 MeV state in 8 Be and not e.g. the ground state 
will be explained below. 

The reaction rate depends critically on the resonance 
energy, E r , i.e. the energy of the 16.8 MeV state relative 
to the d + 7 Be threshold at S d = 16.4901(10) MeV 0]: if 
too far above the threshold, the tunneling process will 
be too slow at BBN temperatures. Furthermore, for 
the proposed destruction mechanism to be efficient, the 
16.8 MeV state must have an appreciable width, Td, for 
being formed in the d + 7 Be channel, but also an appre- 
ciable width, T — Td, for not decaying back to d + 7 Be. 
The energetically allowed decay modes competing with 
deuteron emission are: 7, p, a, and 3 He. However, 7 
and 3 He can safely be neglected. A deuteron width, Td, 
of the required magnitude can only be realized if the 
16.8 MeV state is not too close to the threshold. The 
analysis of Ref. Q shows that the cosmological lithium 
problem can be resolved provided E r ~ 170 — 220 keV, 
T d ~ 10 - 40 keV, and r - T d ~ T d . At the time Ref. 
was written, the known properties of the 16.8 MeV state 
did not contradict these requirements: The 16.8 MeV 
state had been observed in two experiments [HI, HH . Its 
energy had been determined to be 16.7 MeV with an un- 
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certainty of 100 keV, and only an upper limit of 100 keV 
existed on its total width. Its spin and parity had not 
been determined, though a tentative 5/2+ assignment 
had been made [13J based on comparison to the mirror 
nucleus, 9 Be. No information existed on its decay prop- 
erties. 

As noted in Ref. 0], the simultaneous requirement of 
E r ~ 170 - 220 keV and r d ~ 10 - 40 keV is physi- 
cally possible, but implies some rather special properties 
for the 16.8 MeV state: a reduced deuteron width com- 
parable to the Wigner limit and a very large channel 
radius of at least 9 fm. In addition, the proposed de- 
struction mechanism could only be reconciled with the 
direct measurement of Ref. [14( with considerable diffi- 
culty: the proton and a decay of the 16.8 MeV state 
had to be dominated by a single proton-decay branch to 
the 16.626 MeV, 2 + state in 8 Be, because decays to the 
lower-lying states in 8 Be would have produced protons 
of sufficient energy to be detected by the experimental 
setup of Ref. [ij]. 

Slightly above the 16.626 MeV state, at 16.922 MeV, 
lies another 2 + state in 8 Be. The two are nearly maxi- 
mally mixed in isospin (I) [l5| and are often referred to 
as the 2^ doublet. The structure of the 16.626 MeV state 
is primarily that of a \p\i<i proton orbiting a 7 Li core in 
its ground state, and the structure of the 16.922 MeV 
state is primarily that of a lpi/2 neutron orbiting a 7 Be 
core in its ground state [l6j]. The analysis in Ref. [l3j 
suggests that the structure of the 16.8 MeV state in 9 B 
is primarily that of a 2s\/2 proton orbiting the 1=1 
component of the 2 + doublet. As noted in Ref. Q, this 
provided reason to think that the overlap between the 
16.8 MeV state in 9 B and p + 8 Be* might indeed be con- 
siderably larger for the 16.626 MeV state than for any of 
the lower- lying states in 8 Be. 

The nuclear physics input for BBN calculations lias re- 
cently been reexamined in Ref. (l7j . which includes new 
reactions, studies the potential effects of reactions for 
which data do not exist, studies the effects of non-thermal 
particles (highly energetic particles produced in exother- 
mic reactions), in particular neutrons which take much 
longer to thcrmalize than charged particles, and includes 
thermal excitation of the first excited states in Li and 
7 Bc. Ref. [l7| concludes that there is little chance of solv- 
ing the cosmological lithium problem with conventional 
nuclear physics, but retains the destruction mechanism 
proposed by Ref. Q as an "alluring" possibility. 

A very recent paper [l8[ reports on a new 9 Be( 3 He,i) 9 B 
measurement performed with the purpose of studying 
Gamow- Teller transition strengths in the A = 9 sys- 
tem. The beam energy was 140 MeV/nucleon, and tri- 
tons were detected in a high-resolution spectrometer at 
scattering angles around 0°. The excitation energy res- 
olution achieved was 30 keV. The 16.8 MeV state is 
strongly excited, and its energy and width are determined 
to be 16.800(10) MeV and 81(5) keV, respectively, in 
good agreement with the two previous experiments. The 
nearby J 7 * = l/2~, I = 3/2 state at 17.076(4) MeV is 



also strongly excited and its energy is well-known from its 
7 decay to the ground state. This gives strong confidence 
in the new energy determination. It is also worth noting 
that the observed angular distribution of the 16.8 MeV 
state is consistent with the proposed 5/2 + assignment. 

To assess the consequence of the new experimental 
data for the proposed destruction mechanism, we em- 
ploy the standard Kawano/ Wagoner BBN code (l9l. [20|. 
We modify the 7 Be(d,p) reaction rate by adding the extra 
term HH 

N A (av) = iV A (8/7r/i 27 ) 1/2 (feT)- 3 / 2 

x / E a(E) exp (—E/kT) dE , (1) 
Jo 

where Na is Avogadro's constant, fi 2 7 is the reduced 
mass, k is Boltzmann's constant, T is the temperature, 
E is the relative kinetic energy, and cr(E) is the cross sec- 
tion for d + Be — > 9 B* —> p + 8 Be*, given by the single- 
channel, single-level approximation of the i?-matrix the- 
ory HI]: 



a{E) = ttA 2 
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where A = h/p = h/ [2^,27 E) 1 / 2 , and ui is a statistical 
weight factor that depends on the spins involved, 
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where J = 5/2 is the (assumed) spin of the 16.8 MeV 
state, ji = 1 is the spin of the deuteron, and j2 = 3/2 
is the spin of 7 Be. Furthermore, T = T 7 + T p + Yd + 
r3 Ho + r a is the total width. We assume T 7 , r3 Hc , and 
r Q to be negligible and T p to be dominated by the decay 
to the 16.626 MeV state in 8 Be. As the 16.8 MeV state 
is located close to threshold, the energy dependence of 
the deuteron width must be taken into account [23| : 

T d = 2P l=l {E) 1 2 d . 

Similarly, for the proton width: 

T p = 2 P t=0 {E') j 2 p , (2) 

where Pi is the penetrability, £ is the orbital angular mo- 
mentum, 7d (jp) is the deuteron (proton) reduced width, 
and E' is the p + 8 Be* relative kinetic energy, 

E' = E + S d - S p - 16.626 MeV 

with S p = -0.1851(10) MeV We note that Eq. © 
is only approximately valid, as it assumes that the 
width of the 16.626 MeV state in 8 Be can be neglected, 
whereas the state actually has a considerable width of 
108.1(5) keV Q with an asymmetric line shape owing to 
interference with the 16.922 MeV state. Still, the approx- 
imation is adequate for the present analysis. Finally, the 
shift, A, is calculated as 

A = -(S e=1 (E)-B)jj ~(S e = (E')-B')j 2 , 
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FIG. 2: Maximal resonant contribution of the 16.8 MeV state 
in 9 B to the 7 Be(d,p) reaction rate calculated from Eq. |T|) 
versus temperature. For comparison, the rate calculated using 
the narrow resonance approximation is also shown. 



where St is the shift function, and the boundary condi- 
tions are B = S l=1 {E r ) and B' = St=n(E' r ). The defi- 
nitions of Pi and Si are given in Ref. [22(. To evaluate 
Pi and Se, suitable channel radii, 027 and ais, must be 
chosen for the formation and destruction channel. 

Relying on the data from the new 9 Bc( 3 Hc,£) 9 B mea- 
surement [lH, we use E r = 310(10) kcV for the reso- 
nance energy and T° = 81(5) keV for the total width. 
The superscript refers to the value at resonance en- 
ergy, i.e. at E = E r . To maximize the reaction rate, we 
chose 7^ and 7 P such that Y° d = T° = 0.5 T . We do 
not have complete liberty in our choice of 7^ and 7 P as 
they should not exceed the corresponding Wigner lim- 
it^ 7w,rf = 3ft 2 /( 2 M27a| 7 ) and 7w, P = 3h 2 /{2^ 18 al 8 ). 
We find that a permissible choice of jd only exists for 
a27 > 6.5 fm, whereas a permissible choice of 7 p exists 
for essentially any value of aig. For the present calcu- 
lation, we chose 027 = 7 fm and a±g = 5 fm, yield- 
ing 7d/7w,d = °- 93 and 7p/7w, P = °- 09 - As argued 



in Ref. 0, 027 = 7 fm represents a physically plausible 
choice of channel radius. The temperature dependence 
of the reaction rate is shown in Fig. [2] For comparison, 
we also show the rate obtained in the narrow resonance 
approximation. We find that the reduction in 7 Li abun- 
dance caused by the inclusion of the resonant contribu- 
tion of the 16.8 MeV state in 9 B to the 7 Be(d,p) reac- 
tion rate is at most 3.5(8)%. This result is essentially 
independent of the choice of channel radii. The quoted 
uncertainty mainly reflects the 10 keV uncertainty on the 
energy determination of the 16.8 MeV state with a small 
contribution (0.2%) from the 5 keV uncertainty on the 
width determination. 

We stress that the assumption of a dominant proton- 
decay branch to the 16.626 MeV state is, by no means, 
important to the conclusion of the present analysis. If 
the decay is assumed to proceed by proton emission to 
lower-lying states in 8 Be or a emission to 5 Li, a similar 
reduction in 7 Li abundance is obtained. The assumption 
of a dominant proton-decay branch to the 16.626 MeV 
state was made mainly to avoid conflict with the direct 
measurement of Ref. [Tj] . 

In summary, we have shown that the 16.8 MeV 
state in 9 B is unable to enhance the 7 Be(d,p) reaction 
rate by the amount needed to resolve the cosmological 
lithium problem. With the new precise determination 
of the energy of the 16.8 MeV state [l8|, the reduc- 
tion in 7 Li abundance owing to the inclusion of the 
resonant contribution of the 16.8 MeV state to the 
7 Bc(d,p) reaction rate is at most 3.5(8)% and probably 
much lower depending on the decay properties of the 
16.8 MeV state which remain unknown. In line with 
Ref. [I?}, we conclude that all possibilities for solv- 
ing the cosmological lithium problem by conventional 
nuclear physics means now seem to have been exhausted. 
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